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Hypervelocity stars (HVS) travel with velocities so high, that they exceed the 
escape velocity of the Galaxy ( 1 - 3 ). Several acceleration mechanisms have been 
discussed. Only one HVS (US 708, HVS 2) is a compact helium star ( 2 ). Here 
we present a spectroscopic and kinematic analysis of US 708. Travelling with 
a velocity of ~ 1200kms _1 , it is the fastest unbound star in our Galaxy. In 
reconstructing its trajectory, the Galactic center becomes very unlikely as an 
origin, which is hardly consistent with the most favored ejection mechanism 
for the other HVS. Furthermore, we discovered US 708 to be a fast rotator. 
According to our binary evolution model it was spun-up by tidal interaction in 
a close binary and is likely to be the ejected donor remnant of a thermonuclear 
supernova. 

According to the widely-accepted theory for the acceleration of HVS, a close binary is dis¬ 
rupted by the supermassive black hole (SMBH) in the centre of our Galaxy and one component 
is ejected as hypervelocity star ( 4 ). In an alternative scenario US 708 was proposed to be ejected 
from an ultra-compact binary star by a thermonuclear supernova type la (SNIa) ( 5 ). However, 
previous observational evidence was insufficient to put firm constraints on its past evolution. 
Here we show that US 708 is the fastest unbound star in our Galaxy, provide evidence for the 
supernova ejection scenario and identify a progenitor population of SN la. 

In contrast to all other known HVS US 708 has been classified as hot subdwarf star (sdO/B). 
Those stars are evolved, core helium-burning objects with low masses around 0.5 M 0 . About 
half of the sdB stars reside in close binaries with periods ranging from ~ 0.1 d to ~ 30 d 
( 6 , 7 ). The hot subdwarf is regarded as the core of a former red giant star that has been stripped 
off almost all of its hydrogen envelope through interaction with a close companion star ( 8 , 9 ). 
Apparently single hot subdwarf stars like US 708 itself are known as well. However, binary 
evolution has also been proposed in this case, since the merger of two helium white dwarfs 
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(He-WDs) is a possible formation channel for those objects ( 10 ). 

The hot subdwarf nature of US 708 poses a particular challenge for theories, that aim at 
explaining the acceleration of hypervelocity stars. Within the slingshot scenario proposed by 
Hills, a binary consisting of two main sequence stars is disrupted by the close encounter with 
the SMBH in the centre of our Galaxy. While one of the components remains in a bound 
orbit around the black hole, the other one is ejected with high velocity ( 4 ). This scenario not 
only explains the existence of the so-called S-stars orbiting the supermassive black hole in the 
Galactic centre and providing the most convincing evidence for the existence of this black hole 
( 11 ). It is also consistent with the main properties of the known HVS population consisting 
of young main sequence stars ( 12 , 13 ). However, more detailed analyses of some young HVS 
challenge the Galactic centre origin ( 3 , 14 ) and most recently, a new population of old main 
sequence stars likely to be HVS has been discovered. Most of those objects are also unlikely to 
originate from the Galactic centre, but the acceleration mechanism remains unclear ( 15 ). 

In the case of the helium-rich sdO (He-sdO) US 708 the situation is even more complicated. 
In contrast to all other known HVS, which are normal main sequence stars of different ages, 
this star is in the phase of shell helium-burning, which only lasts for a few tens of Myr. More 
importantly, it has been formed by close binary interaction. To accelerate a close binary star to 
such high velocity with the slingshot mechanism, either a binary black hole ( 16 ) or the close 
encounter of a hierarchical triple system, where the distant component becomes bound to the 
black hole and the two close components are ejected, is necessary ( 17 ). Similar constraints ap¬ 
ply to the dynamical ejection out of a dense cluster, which is the second main scenario discussed 
to explain the HVS. 

While close binarity requires quite specific modifications of the canonical HVS scenarios, 
it is a necessary ingredient for an alternative scenario, where US 708 is explained as the ejected 
donor remnant of a thermonuclear supernova type la (SNIa) ( 18 , 19 ). Underluminous SNIa 
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have been proposed to originate from a so-called double-detonation of a white dwarf ( 20 , 21 ). 
In this scenario a massive white dwarf (WD) is closely orbited by a low-mass helium star. Due 
to a tightening of the orbit the He-star will start to transfer mass to its compact companion. 
After a critical amount of helium is deposited on the surface of the WD through accretion, the 
helium is ignited causing a detonation wave that triggers the explosion of the carbon-oxygen 
white dwarf itself. Indeed, the ultracompact sdB+WD binary CD—30°11223 has recently been 
identified as progenitor candidate for such a scenario and linked to the putative ejected donor 
remnant US 708 ( 5 , 22 ). 

We performed a detailed spectroscopic and kinematic analysis of US 708 based on recently 
obtained and archival data to trace back its origin and constrain the ejection mechanism. To 
determine the 3D-motion of US 708 both the radial and tangential velocity components have 
to be determined. We measured the radial velocity from new spectra taken with the Keck and 
Palomar telescopes and compared it to archival data. Atmospheric parameters and spectro¬ 
scopic distance were derived from the new spectra as well (see Fig. 1). In addition to that the 
proper motion was determined by combining archival positions with new measurements from 
the PanSTARRS survey (see Fig. S2, Table 1). 

With a Galactic restframe velocity of 1157 ± 53kms _1 we find that US 708 is the fastest 
known unbound star in our Galaxy. Its current distance is 8.5 ± 1.0 kpc and it is moving away 
from the Galactic plane into the halo. Tracing back its trajectory and assuming no further 
deviations we deduced that it crossed the Galactic disc 14.0 ± 3.1 Myr ago. In this way an 
origin in the centre of our Galaxy can be excluded with high confidence (see Fig. 2), but the 
origin in the Galactic disc on the other hand is fully consistent with the supernova ejection 
scenario. In contrast to regular SNIa, double-detonation SNIa with hot subdwarf donors are 
predicted to happen in young stellar populations ( 5 ). 

Both the current Galactic restframe velocity and the reconstructed ejection velocity from the 
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Table 1: Parameters of US 708. The uncertainty of the radial velocity is the la error from 
a y 2 -fit. the uncertainties in the proper motion components have been propagated from the 
position errors by linear regression and the uncertainties in the atmospheric parameters are 
bootstrap errors. The uncertainties of the other, derived parameters have been propagated from 
the uncertainties of the input parameters. 


Visual magnitude^ 

m g 

[mag] 

18.668 ±0.008 

Proper motion 

/la COS (5 

[mas/yr] 

-8.0 ±1.8 


Vs 

[mas/yr] 

9.1 ± 1.6 

Radial velocity 

Gelio 

[kms -1 ] 

917 ±7 

Galactic cartesian coordinates 

X 

[kpc] 

-14.2 ±0.7 


Y 

[kpc] 

0.4 ±0.1 


Z 

[kpc] 

6.2 ±0.7 

Galactic velocity components 

vx 

[kms -1 ] 

-855 ±61 


Vy 

[kms -1 ] 

643 ± 77 


vz 

[kms -1 ] 

431 ± 58 

Galactic rest-frame velocity 

V gr f 

[kms -1 ] 

1157 ±53 

Effective temperature 

T e fi 

[K] 

47 200 ±400 

Surface gravity 

l°g g 


5.69 ±0.09 

Helium abundance 

log y 


±2.0 

Nitrogen abundance 

be im 

iU & JV(H) 


-2.4 ±0.2 

Projected rotational velocity 

Got sin i 

[kms -1 ] 

115 ±8 

Mass (adopted) 

M S dB 

\M & \ 

0.3 

Distance 

d 

[kpc] 

8.5 ±1.0 


t Taken from the Sloan Digital Sky Survey Data Server (das.sdss.org) 


5 



Galactic disc (998 ± 68kms -1 ) are significantly higher than published before (~ 750 km s -1 , 
based on radial velocity alone) (2). This puts new constraints on the possible progenitor system, 
which can be derived from the observed parameters of US 708. To reach such a high ejection ve¬ 
locity, the progenitor binary must have been very compact and the white dwarf companion rather 
massive. The likely progenitor system consisted of a compact helium star with a mass of ~ 
0.3 M 0 and a massive carbon-oxygen WD (1.0 —1.2 M 0 ) with an orbital period of about 10 min. 
We calculated the mass-transfer rate in such a binary and found that the helium is accreted by 
the white dwarf at a rate suitable for the double-detonation scenario (10' 9 — 10" 8 M 0 yr _1 , for 
details see (5)). Such ultra-short period systems with compact helium stars have indeed been 
observed. The eclipsing He-WD+CO-WD binary SDSS J065133+284423 has an orbital period 
of only 12 min (23). However, the mass of the CO-WD (0.55 M 0 ) is too low for a double¬ 
detonation SN la. 

The ejection from such a close binary should leave another imprint on the remnant. We know 
that hot subdwarfs in compact binaries have been spun up by the tidal influence of the close 
companion (22,24,25) to rotational velocities significantly higher than the rotational velocities 
of single hot subdwarfs (26,27). An ejected remnant is predicted to have a high u rot sin i as well 
(28). And indeed we measured v mt sin i = 115 ± 8kms _1 significantly higher than expected 
for a single He-sdO (see Fig. 1) (27). However, assuming an ultracompact progenitor binary, 
the measured i; ro1 sin i of the He-sdO is still about a factor of four slower than expected. A 
significant loss of mass and angular momentum caused by the supernova impact is predicted 
for main sequence companions in the standard single-degenerate SN la scenario, where mass 
is transferred from such a companion to a WD (29-31). While the loss of mass and angular 
momentum for a more compact companion like US 708 is expected to be smaller, our results 
indicate that it is still substantial. This puts important observational constraints on more detailed 
models. 
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While the observed properties of US 708 are consistent with the supernova ejection sce¬ 
nario, they are hardly compatible with the slingshot mechanism, because an origin of the star in 
the centre of the Galaxy is very unlikely (Fig. 2, see also the additional explanation in the sup¬ 
plementary material). However, it has to be pointed out that the supernova ejection scenario is 
only applicable to such compact helium stars and cannot be invoked to explain the acceleration 
of the other HVS. 

Providing evidence that US 708 is the likely donor remnant of a helium double-detonation 
SN la we not only show that the fastest unbound stars in our Galaxy are accelerated in this way. 
It is also an important step forward in our understanding of SN la explosions in general. Despite 
the fact that those bright events are used as standard candles to measure the expansion (and 
acceleration) of the universe, their progenitors are still unknown. Our results suggest that due 
to the quite high WD masses derived for the progenitor binary, the double-detonation scenario 
might not only be applicable to some underluminous SNIa (5,21), but might also contribute to 
the population of typical SNIa used as cosmic yardsticks (20). 

Depending on the pollution by SN material, the effect of the SN impact, and the subse¬ 
quent stellar evolution, the surface abundances of US 708 might be significantly affected. UV- 
spectroscopy is necessary to measure the metal abundances of US 708 and put further con¬ 
straints on the extreme history of this star, which witnessed a SN from a distance of less than 
0.2 R e . 
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Fig. 1: Fit of model spectrum. Fit of synthetic NLTE models to the helium and nitrogen 
lines of a Keck/ESI spectrum of ETS 708. The normalized fluxes of the single lines are shifted 
for better visualisation and the most prominent lines are labeled. The weaker lines are from 
N ill at 4634 A, 4640 A, He I at 4713 A (middle plot), N til at 4379 A and He I at 4387 A (bottom 
plot). The dashed vertical line marks the rest wavelengths of the lines. The high radial velocity 
shift as well as the significant broadening of the lines are clearly visible. 


9 












I 20000 
I 10000 

! 15000 

H 2000 

t 4 iooo 

1 - 500 

- - 200 
C 100 

; la 


= 

s 


s- 

<v 

a 


a 

<u 

E 

•- 

3 

ui 

S3 

<u 


Cm 

O 

i- 

<v 

s 

3 


X in kpc 

Fig. 2: Origin of US 708. Monte Carlo simulation (10 8 iterations) of the past trajectory of 
US 708. The colour-coded bins mark the positions, where the star crossed the Galactic disc, 
which is shown pole-on. The contours correspond to the la, 3a and 5a confidence limits. The 
position of the Galactic centre is marked by the black dot, the position of the Sun with the star 
symbol. The current position of US 708 is marked by a triangle and given in Table 1. 
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Materials and Methods 

Summary 

Spectra of US 708 have been taken with the 10m Keck and the 5m Palomar telescopes. From the 
Doppler shift of the spectral lines we measured the radial velocity using both new and archival 
data. The tangential velocity components have been determined by measuring the proper motion 
of the star from multi-epoch position measurements spanning 59 years and its spectroscopic dis¬ 
tance performing a full quantitative spectral analysis using state-of-the-art model atmospheres. 
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Using those informations we constrained the kinematics of this star and traced back its origin 
to the Galactic disk performing a Monte Carlo simulation. The properties of the progenitor 
binary have been determined mostly based on the derived ejection velocity from the Galactic 
disc. Binary evolution calculations have then been performed to check the consistency of those 
properties with theory. The current rotational properties of US 708 have been compared with 
hydrodynamical models of angular momentum-loss triggered by supernovae explosions. 

Observations 

US 708 (a 2 ooo = 09 h 33 m 20?85, <5 20 oo = +44°17'05.8") was discovered to be a hypervelocity 
star by Hirsch et al. (2). A medium-resolution ( R ~ 1800) spectrum was taken by the Sloan 
Digital Sky Survey (SDSS) on February 20, 2002 (32). Follow-up low-resolution ( R ~ 900) 
spectroscopy was obtained with the LRIS instrument at the Keck telescope on May 13, 2005. 
The reduced spectra from the blue and red channel of the instrument were provided to us by 
H. Hirsch. A series of 11 consecutively taken medium-resolution ( R ~ 8000) spectra was 
obtained with the ESI instrument at the Keck telescope on March 3, 2013. The spectra have 
been reduced with the ESI pipeline Makeefl One spectrum has been taken with the medium- 
resolution spectrograph at the 5.1 m Hale telescope on Mount Palomar on May 11, 2013 and 
another three spectra on June 1, 2013. 

Revised radial velocity 

Hirsch et al. measured the radial velocity (RV) of US 708 (708 ± 15 kms -1 ) from the helium 
lines in the blue-channel LRIS spectrum. The measured RV exceeded the typical RVs of He- 
sdOs in the rest of the sample, which is characteristic for halo stars, considerably (see Fig. 5 in 
(35)). We obtained the RV of the SDSS, ESI and Palomar spectra by fitting a model spectrum 

1 http://www.astro.caltech.edu/~tb/makee/ 
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(see below) to the helium lines using the FITSB2 routine (33). Surprisingly, the most accurate 
RV measured from the coadded ESI spectrum (917 ± 7 kms -1 ) turned out to be significantly 
higher than the one published by Hirsch et al. (2). This velocity is consistent with the RVs 
measured both from the SDSS^I (898 ± 30 km s -1 ) and the Palomar spectra (866 — 936 km s -1 ). 
To investigate this issue, we performed a reanalysis of the LRIS spectra and measured an RV 
of 709 ± 7 kms -1 for the LRIS blue-channel spectrum perfectly consistent with the published 
value. However, when fitting the red-channel spectrum we found a significantly discrepant RV 
of 797 ± 21 kms -1 . This was taken as first indication, that those spectra might be affected 
by systematics. We used the nightsky emission line of Ol at 6300 A (red-channel) and the 
interstellar absorption line of Ca II at 3934 A (blue-channel) to quantify those systematic shifts 
(see Fig. SI, left panel). The nightsky emission line, which is supposed to be at zero RV, 
was blue-shifted by —33 ± 10kms -1 . The interstellar line showed a significantly higher shift 
of —128 ± 22kms -1 . Since the correct RV of the interstellar line is not known a priori, we 
measured it from a coadded Palomar spectrum to be —9 ± 28kms -1 . Correcting the RVs for 
those shifts, the two RV values (828 ± 22kms -1 blue-channel, 830 ± 21 kms -1 red-channel) 
from the LRIS spectra are consistent with each other, but still smaller than the RVs measured 
from the ESI and Palomar spectra (see Fig. SI, right panel). Due to the low resolution of 
the LRIS spectra, the remaining shift corresponds to only about one pixel on the CCD and is 
therefore regarded as systematic as well. We conclude that the RV published by Hirsch et al. 
was affected by systematics and therefore underestimated. Going back to the original raw data, 
we performed an independent reduction. However, we were not able to resolve this issue. Given 
that the uncertainties are at the lcr level of confidence and that the LRIS spectra are affected by 
systematics, no significant shifts in RV on both short and long timescales are detected (see 

2 The RV of 793 kms -1 given by the SDSS Sky Server Object Explorer tool is based on the fit of an hydrogen 
rich template to the spectrum. The He Il-lines of the Pickering series are misidentified as Balmer lines. This 
introduces the shift of ~ —100 km s -1 between our result and the template fit. 
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Supplementary Fig. 1, right panel). 


Proper motion 


The proper motion of US 708 has been derived from multi-epoch position measurements of 


Schmidt plates obtained from the Digitised Sky Survey (DSS 


3 , the Sloan Digital Sky Survey 


(32) and the PanSTARRS survey (PS1) over a timebase of 59 years (see Fig. S2). The posi¬ 
tions from the DSS and SDSS images have been measured with respect to a set of compact 
background galaxies selected from the SDSS images as described in Tillich et al. (34). The 
positions of the background galaxies and the object are measured. For each image the measured 
positions of the background galaxies are compared to the reference values from the PS 1 cata¬ 
logue. The average of the deviations from these reference values is adopted as uncertainty of the 
object position. In the case of the 29 PS1 epochs, we took the calibrated positions from the PS1 
catalogue and therefore use the same reference system for all our measurements. We obtained 
one position per epoch and used linear regression ro derive the proper motion components with 
their uncertainties. 


Atmospheric parameters 

The atmospheric parameters effective temperature T eff , surface gravity log g, nitrogen abun¬ 
dance log N(N)/N(H.) and projected rotational velocity ?; rot si i 1 1 were determined by fitting 
simultaneously the observed helium and nitrogen lines of an ESI spectrum, constructed by 
coadding the 11 single exposures, with NLTE models taking into account line-blanketing of ni¬ 
trogen (27) (see Fig. 1) as described in Geier et al. (35). Since no hydrogen lines are visible in 
the spectrum, we fixed the helium abundance to log y = log A^(He)/iV(H) = +2.0. The atmo¬ 
spheric parameters (T eff = 47200 ± 400 K, log g = 5.69 ± 0.09) deviate significantly from the 

; http ://arch ive.stsci.edu/cgi-bin/dss_plate_fin tier 
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results by Hirsch et al. (T e g = 45600 ±700 K, log g = 5.23±0.12) especially in surface gravity. 
This is caused by the additional line-blanketing of nitrogen (27). The atmospheric parameters as 
well as the nitrogen abundance log 7V(N)/1V(H) = —2.4 ± 0.2 are typical for the nitrogen-rich 
subclass of He-sdOs (36). The projected rotational velocity v Iot sin i = 115 ± 8 km s _1 on the 
other hand is significantly higher than the ones of both single sdB (< lOkms -1 ) and He-sdO 
stars (20 — 30kms _1 ) (26,27). Based on our analysis we can rule out objects with similar 
spectral features like DO-type white dwarfs (log g > 7.0) or luminous He-stars (log g < 4.5), 
which can be easily misclassified from visual inspection only. 


Spectroscopic distance and kinematics 


The spectroscopic distance is derived from the atmospheric parameters T eS , log g and the appar¬ 
ent visual magnitude as described in Ramspeck et al. (40). The SDSS-g and r magnitudes have 
been converted to Johnson V magnitude!, which has been corrected for interstellar reddening 
(A\i = 0.07 mag) (41). Based on constraints provided by the supernova ejection scenario (see 
below) we adopted a mass of 0.3 M 0 for the hot subdwarf. The spectroscopic distance in this 
case is 8.5 ± 1.0 kpc. For this distance the proper motion components are converted to absolute 
transverse velocities and, combined with the radial velocity, the Galactic restframe velocity of 
US 708 is calculated to be 1157 ± 53kms _1 . This is the highest known restframe velocity of 
any unbound star in our Galaxy. The past trajectory of US 708 in the Galactic potential (42,43) 
has been reconstructed as outlined in Tillich et al. (34). Due to the high velocity of US 708 
we found that the trajectory is not changed significantly if alternative model potentials are used 
(42). US 708 was ejected from the Galactic disc 14.0 ±3.1 Myr ago and the ejection velocity, 
corrected for the motion of the Galactic disc, was calculated to be 998 ± 68kms _1 . We per¬ 
formed Monte Carlo simulations to trace back the trajectory of US 708 until it intersects with 


4 http://www.sdss.org/dr6/algorithms/sdssUB VRITransform.html 
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the Galactic disc. The uncertainty in the proper motion measurement dominates the error bud¬ 
get. Assuming no further perturbation of the trajectory an origin from the central kpc around 
the Galactic centre can be ruled out with a confidence of more than 6cr (see Fig. 2). 

Properties of the progenitor binary. Geier et al. (5) proposed the ultra-compact sdB+WD 
binary CD—30° 11223 to be a possible progenitor of the hypervelocity sdO US 708. However, 
based on the new results presented here, the ejection velocity is ~ 250kms _1 higher than 
assumed by Geier et al. Hence, it is necessary to reexamine the supernova ejection scenario 
and to test its consistency with the newly derived parameters of US 708. Similar to the scenario 
discussed in Geier et al., we assume that the progenitor binary consisted of a compact helium 
star and a massive carbon-oxygen WD in close orbit. The ejection velocity of the He-star 
equals the radial velocity semiamplitude of the progenitor binary at the moment of the supernova 
explosion (K = 998 ± 68 km s _1 ) modified by the additional perpendicular velocity component 
the star received through the SN explosion (~ 200 km s _1 ). Since both velocities are added in 
quadrature, the kick velocity is negligible within the uncertainties. Assuming a mass for the 
He-star and a circular orbit, the orbital period of the progenitor binary as well as its separation 
can be calculated from the binary mass function: 


fm = 


^WD 


sin 3 i 


PK 3 


(M W d + Mn e ) 2 2nG 

Since we known the absolute space velocity, the inclination angle can be set to sin i 


( 1 ) 

= 1 , 


and the orbital period P of the progenitor binary can be calculated: 


P 


2v tG 


^WD 


(M W d + M] 


HcJ 


The binary separation a can be derived using Keplers laws: 


( 2 ) 
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a = 


G M~wd 


(3) 


A 2 M wd + M He 

Another crucial assumption is that stable mass-transfer from the He-star to the WD trig¬ 
gered the SN, which means that the He-star must have filled its Roche lobe before ejection. 
To calculate the Roche lobe radius we used the equation given by Eggleton (44), where q = 

M He /M MwD : 


At = 


0.49g 2 / 3 c 


(4) 


0.6g 2 / 3 + ln(l + q 1 / 2 ) 

The radius of the He-star can be calculated as a function of the mass and the surface gravity 


9- 


AHe ~ 



(5) 


To compare the Roche radius with the possible radius of the He-star we have to take into 
account that US 708 has already evolved away from the EHB (see Fig. S3), which led to an 
increase in radius. To calculate the radius at the time of ejection we therefore adopt the highest 
reasonable values for log g ~ 6.1 close to the ZAEHB and the He-MS. Calculating Roche lobe 
and He-star radii for different He-star and WD masses, we explored the parameter space and put 
constraints on possible progenitor systems. Fig. S4 shows the Roche radii for WD masses from 
1.0 M e to 1.2 M q . The He-star radii for log g = 6.1 are plotted for comparison. Consistent solu¬ 
tions are only found for low He-star masses (~ 0.3 — 0.35 M & ) and high WD companion masses 
(~ 1.0 — 1.2 M q ). The orbital period of the progenitor binary can be constrained to ~ 10 min. 
Following the method described in Geier et al. (5) we calculated the mass-transfer rates for 
binaries with similar orbital parameters and component masses. The rates of 10“ 8 M 0 yr" 1 
are consistent with the helium double-detonation scenario. Figs. S5-S7 show as an example the 
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evolution of a close binary (P = 26 min) that starts mass-transfer with an initial He-star mass 
of 0.45 Mq and a WD mass of 1.05 M 0 . After about 5 Myr the orbital period becomes as short 
as 13 min and the component masses change to 0.3 M 0 and 1.2 M & , when the WD explodes as 
SN la. The helium donor has to be a helium-burning star rather than a He-WD without ongoing 
nuclear burning in its core, because those objects consist of degenerate matter and as soon as 
the small non-degenerate envelope is transferred, the mass-transfer rate becomes too high for 
the double-detonation scenario. Such systems will experience He-flashes on the surface of the 
WD companion without igniting the carbon in the core (45). This is consistent with the obser¬ 
vational evidence. Since the ejection already happened ~ 14 Myr ago, we can also assume that 
US 708 is a helium-burning star. The minimum mass for such objects is ~ 0.3 Mq. Even less 
massive He-stars without nuclear burning, which are the progenitors of He-WDs, exist (46). 
However, according to evolutionary tracks, their effective temperatures are much lower than 
the one of US 708. The most massive He-WD progenitors on the other hand, which can reach 
such temperatures, are cooling within a few Myr, too fast to be consistent with the position 
of US 708 in the T e g — log ^-diagram (see Fig. S3) (47). EHB-tracks for masses as low as 
~ 0.3 M q are also not consistent with the position of US 708 in the T eS — log ^-diagram (6). 
However, since the He-star was significantly more massive before the mass-transfer started, its 
further evolution might not depend on its current total mass. Especially, if the helium in its core 
was already exhausted towards the end of the mass-transfer phase, the further evolution would 
depend on the core mass rather than the total mass. The position of US 708 in the T eS — log g- 
diagram is, for example, perfectly consistent with post-EHB model tracks for an original mass 
of 0.45 M & (see Fig. S3). Based on these simple calculations we can rule out the sdB+WD 
binary CD—30° 11223 (0.51 M 0 + 0.76 Mq, P = 72 min) as direct progenitor of US 708. How¬ 
ever, systems like CD—30° 11223 remain progenitor candidates of other high velocity sdB stars 
(34). There is no binary known yet, which fulfills all the criteria for a progenitor system to 
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US 708. But evidence is growing, that such objects exist. A whole population of close binaries 
consisting of He-WDs and CO-WDs has been discovered recently (46). They form in almost 
exactly the same way as sdB+WD binaries. The only difference is that core helium-burning was 
not ignited before the hydrogen envelope has been stripped off in the common envelope phase 
of unstable mass-transfer. The ultracompact, eclipsing He-WD+CO-WD (0.25 M 0 + 0.55 M 0 ) 
binary SDSS J065133+284423 with an orbital period of only 12 min sticks out (23). Its period 
is similar to the one expected for the progenitor of US 708, but the component masses are too 
small. In the double-lined WD+WD binary SDSS J125733+542850 (0.2 M 0 + 1.2 M 0 ) on the 
other hand, the masses are very similar to the ones predicted, whereas the orbital period is much 
longer (274 min) (48). These discoveries as well as binary evolution calculations indicate the 
existence of binaries fulfilling the criteria for a progenitor of US 708 as well (18). We there¬ 
fore conclude that the double-detonation supernova ejection scenario is still able to explain the 
observed properties of US 708 as ejected donor remnant. 

Rotational velocity 

Only two out of more than 100 single hot subdwarf stars are fast rotators. Both objects are 
sdB stars with hydrogen-rich atmospheres and might have been formed by a common-envelope 
merger (49,50). US 708 on the other hand belongs to the population of He-sdOs, which are 
hotter and show no or only some hydrogen in their atmospheres. They are regarded as a distinct 
group of stars, that might have been formed in different ways as the sdBs. US 708 is the only 
single He-sdO rotating faster than 20 — 30 km s -1 indicating a close-binary origin (27). Due to 
the short orbital period and high companion mass, the rotation of the He-star in the proposed 
progenitor binary is expected to be synchronised to its orbital motion (22,24,25). Assuming the 
angular momentum is unchanged after the SN, the ejected donor remnant should remain a fast 
rotator. The rotational velocity can be calculated as follows: 
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( 6 ) 


2nR h 

V IOt — p 

The expected rotational velocity of the ejected He-star is only weakly dependent on the 
masses of the binary components and of the order of 600 km s 1 . This is much higher than the 
measured v TOi sin i = 115±8kms _1 . The significant difference between the expected rotational 
velocity and the measured v TOt sin i comes unexpected. In a synchronised binary system the 
rotational axes of both components are perpendicular to the orbital plane. As soon as the He- 
star is ejected, the rotation axis should be perpendicular to the flight trajectory, which means that 
sin i ~ 1. The impact of the supernova shockwave on main sequence (MS) companions in the 
standard single-degenerate scenario has recently been studied with hydrodynamic simulations. 
Due to stripping of matter, the star may lose up to ~ 90% of its initial angular momentum. A 
subsequent increase in radius due to stellar evolution is also predicted to lower the rotational 
velocity at the surface ( 28 - 31 ). However, simulations of more compact helium stars show that 
much less mass is stripped ( 28 , 51 ). The loss of angular momentum is also expected to be smaller 
in this case. Taking into account evolution on the extreme horizontal branch (EHB), the radius 
of the sdO increased by about a factor of ~ 1.6 since the ejection. Assuming conservation of 
angular momentum, the rotational velocity should now be of the order of 400 km s' 1 . Whether 
the rest of the angular momentum was lost in the SN impact or later is still unclear. Pan et 
al. ( 28 ) predict an increase of the helium star’s radius by a factor of up to four right after the 
impact. This phase should last for a few tens of years. Due to the high initial rotational velocity 
of the star, another episode of mass and angular momentum loss may be possible in this phase. 
Another possibility might be a tilting of the stars rotation axis by the SN impact. The projected 
rotational velocity of the star measured from the line broadening would then be smaller than the 
true rotational velocity. However, simulations show that this effect is negligible for MS stars 
and most likely also for the more compact He stars studied here. 
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Supplementary Text 

Discussing alternative acceleration scenarios 

Any scenario for the acceleration of US 708 must explain four key properties of this star si¬ 
multaneously: (I) US 708 is a compact He-sdO, which most likely formed via close binary 
interaction. Either it is the stripped core of a red giant or the result of a He-WD merger. (II) 
The star has the highest Galactic restframe velocity (~ 1200 km s' ') ever measured for any un¬ 
bound star in the Galaxy. (Ill) The past trajectory of the star is not consistent with an origin in 
the Galactic centre. (IV) In contrast to all other known single He-sdOs, US 708 has a projected 
rotational velocity exceeding ~ lOOkms -1 . We now want to review other scenarios for the 
acceleration of hypervelocity stars in this way. In the classical runaway scenario a massive star 
in a binary system explodes as core-collapse supernova, while the companion is kicked out of 
the system (52). However, the ejection velocity scales with the binary separation and because 
a massive and hence large star is involved, the binary separation cannot become small enough 
to reach an ejection velocity like the one of US 708. The disruption of a hierarchical triple sys¬ 
tem consisting of a normal star in wide orbit around a close He-WD binary by the SMBH is 
regarded as very unlikely as well, because an origin in the GC where the SMBH is located is 
very unlikely. The subsequent merger of a He-WD binary ejected in this way has been proposed 
as formation scenario for US 708 (2,16). The ejection of a He-WD binary star by a hypothetical 
binary black hole in the GC is very unlikely for the same reason (16). Other formation channels 
for hypervelocity stars invoke dynamical interactions in dense stellar clusters (53). Interactions 
of two close binaries can lead to the ejection of a star with the appropriate velocity. However, 
the binary-binary interaction is not affecting the angular momentum of the ejected star. We can 
therefore assume that the observed v rot sin i — 115 ± 8 km s -1 resembles the rotational velocity 
in the tidally-locked progenitor binary. To reach such a high rotational velocity, the separation 


27 



of this binary is constrained to ~ 1 R Q (24). The interaction probability of two such binaries 
even in a very dense cluster is expected to be extremely small. 

Another idea to explain HVS not originating from the GC is the origin in a nearby, low-mass 
galaxy (15). Since the escape velocities from those smaller galaxies are smaller as well, it is 
easier for stars to escape and travel trough the intracluster medium. Some of those neighbouring 
galaxies also have quite high velocities with respect to our own Galaxy. However, this scenario 
is also unlikely for US 708. Although the star might have lived long enough on the main se¬ 
quence (~ 10 Gyr) to travel all the way from a satellite or small neighbouring galaxy, its current 
state of evolution is quite short compared to its total lifetime (only about 0.1%). Furthermore, 
only about 2% of all main sequence stars undergo an EHB phase. This means that for each 
single HVS sdO coming from the intracluster medium there should be about 50 000 HVS main 
sequence stars travelling through our Galactic halo. However, only a few tens of them have 
been reported so far. While faint, high proper motion objects are still not easy to identify, it is 
very easy to discover stars with high RVs in big survey like SDSS or RAVE. Palladino et al. 
(15) list more exotic mechanisms like interactions in globular clusters, with intermediate mass 
black holes or between galaxies (54). In addition to that, combinations of several scenarios 
are imaginable. If for example a hierarchical triple system would be disrupted by the SMBH 
and one component of the ejected binary would explode as core-collapse SN, the trajectory of 
the surviving companion would not point back to the GC. However, we also regard all those 
scenarios as very unlikely to explain the object presented here. 

Estimating supernova rates 

Another sanity check for our scenario is to provide a rough estimate of the double-detonation 
SN la rates we would expect based on our observations and binary population synthesis models. 
So far US 708 is unique among the known He-sdO stars and this estimate is based on very small 
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number statistics. The star was drawn from a sample of hot subdwarfs selected from SDSS. The 
full sample contains 1369 hot subdwarfs in total, 262 of them or roughly 20% are He-sdOs (55). 
Binary population synthesis models by Han et al. (9) predict a birthrate of 5 x 10~ 2 yr -1 for hot 
subdwarfs in general and therefore 1 x 10 2 yr -1 for He-sdOs. One of the observed He-sdOs 
(US 708) might be an ejected donor remnant (~ 0.4%). This translates into a double-detonation 
SN rate of roughly 4 x 10~ 5 yr _1 . This has to be regarded as lower limit only, because a few 
He-sdOs with smaller RVs, but rather high proper motions might still be hidden in our sample. 

The predicted rates of double-detonation SNIa are around 3 x 10~ 4 yr~ 1 and the observed 
rates of all types of SNIa around 3 x 10 -3 yr -1 (56). Since all those numbers have quite 
significant uncertainties, they are regarded as broadly consistent. The most important point for 
this sanity check is, that our estimates from observations do not deviate from the predicted or 
observed rates by orders of magnitude. 

Yu & Tremaine (57) calculated the ejection rates of hypervelocity stars expected from in¬ 
teractions with the Galactic centre black hole (or a binary black hole in the GC). The rates 
are of the order of ~ 10 _5 yr _1 to ~ 10~ 4 yr _1 . However, these numbers correspond to the 
simplest case, the ejection of single main sequence stars. Since stars as peculiar as US 708 and 
their progenitors are very rare compared to normal main sequence stars, the close encounter and 
ejection rates of such stars have to be orders of magnitude smaller. It is therefore very unlikely 
to find one He-sdO along with the about 20 other hypervelocity stars assuming that they are all 
accelerated in the GC. 
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Fig. S1: Revised radial velocity. Left panel: Interstellar Ca II line of US 708 in the Palomar 
(upper plot) and blue-channel LRIS spectra (middle plot). In contrast to the Palomar spectrum, 
the LRIS spectrum is significantly blue-shifted. Nightsky emission line of Ol in the LRIS 
red-channel spectrum (lower panel). The blue-shift is smaller than in the red-channel spec¬ 
trum. Right panel: Radial velocities of US 708 plotted against Julian date. Upper panel: SDSS 
(triangle), LRIS blue- and red-channel uncorrected (grey circles), LRIS blue- and red-channel 
corrected (black circles), ESI (diamonds), Palomar (squares). The dotted box marks the LRIS 
RVs, which are affected by systematics. Lower panel: Close-up of the ESI RVs taken within 
one night. 
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/j a cos6 = —8.0 ± 1.8 mas yr -1 
fig = 9.1 ± 1.6 mas yr -1 



Fig. S2: Proper motion of US 708. Relative positions of US 708 in right ascension (upper 
panel) and declination (lower panel) plotted against time. The POSS I, QUICK and POSSII 
positions are measured from scanned photographic plates provided by the Digitised Sky Survey. 
SDSS and PanSTARRS positions are measured from CCD images. The black solid lines mark 
the best fits, from which we derive the proper motion components. The solid red lines mark the 
proper motion components required for the star to originate from the Galactic centre. 
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Fig. S3: Evolutionary status of US 708. T eS — log g diagram. Evolutionary tracks (solar 
metallicity) of core helium-burning stars with a mass of 0.45 M 0 and different masses of their 
hydrogen envelopes (for bottom to top, 0.0 M 0 , 0.001 M 0 , 0.005 M 0 ) are plotted (6). The 
positions of both the Zero Age and the Terminal Age Extended Horizontal Branch (ZAEHB, 
TAEHB) are indicated as well as the helium main sequence (He-MS). The filled black symbols 
mark known He-sdBs (37,38) (squares) and He-sdOs(39) (triangles) from the literature. 
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Fig. S4: Mass-radius relation of the compact He-star. The dotted black lines mark the Roche 
radii for WD companion masses of 1.0 M e and 1.2 M 0 . The red solid line marks the He-star 
radius assuming log g = 6.1. 
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Fig. S5: Mass-transfer rate. The solid and dash-dotted curves show the mass-transfer rate 
and the mass of the WD envelope (He shell) varying with time after the He star fills its Roche 
lobe, respectively. The dotted vertical line indicates the position where the double-detonation 
may happen (the mass of the He shell increases to ~ 0.15 M 0 ). The initial binary parameters 
and the parameters at the moment of the SN explosion are also given. 
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Fig. S6: Change of component masses. Change of He-star (solid line) and WD mass 
(dash-dotted line) with time. The dotted vertical line indicates the position where the double¬ 
detonation may happen (the mass of the He shell increases to ~ 0.15 M 0 ). The initial binary 
parameters and the parameters at the moment of the SN explosion are also given. 
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Fig. S7: Evolution of orbital parameters. Time evolution of the radial velocity semiampli¬ 
tude (solid line) and the orbital period (dash-dotted line) of the binary. The dotted vertical line 
indicates the position where the double-detonation may happen (the mass of the He shell in¬ 
creases to ~ 0.15 Mq). The initial binary parameters and the parameters at the moment of the 
SN explosion are also given. 
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